Time delays across saddles as a test of modified gravity 
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Modified gravity theories can produce strong signals in the vicinity of the saddles of the total 
gravitational potential. In a sub-class of these models this translates into diverging time-delays for 
echoes crossing the saddles. Such models arise from the possibility that gravity might be infrared 
divergent or confined, and if suitably designed they are very difficult to rule out. We show that 
Lunar Laser Ranging during an eclipse could probe the time-delay effect within meters of the saddle, 
thereby proving or excluding these models. Very Large Baseline Interferometry, instead, could target 
delays across the Jupiter-Sun saddle. Such experiments would shed light on the infrared behaviour 
of gravity and examine the puzzling possibility that there might be well-hidden regions of strong 
gravity and even singularities inside the solar system. 



Even though it is inevitable that Einstein's theory of 
general relativity (GR) will not be the final word, it is 
telling that almost a century after its proposal all the- 
ories trying to supersede it have been ruled out or re- 
main beyond detection Nonetheless, it is precisely 
the experimental misfortunes of "modified gravity" that 
prove the strength of GR, so it is important to keep push- 
ing the boundaries, constructing and observationally dis- 
proving new possibilities. A further motivation derives 
from attempts to combine quantum theory and general 
relativity, a logical (if not an empirical) necessity. Such 
efforts invariably lead to corrections to GR, often at en- 
ergy scales beyond the reach of current experiment, but 
not always. Finally, we should never forget that face 
value — taking into account only the matter sources that 
we do see — the observational status of GR in astrophysics 
and cosmology is calamitous. This is usually blamed on 
our imperfect knowledge of the matter content of the 
Universe, and dismissed by introducing new forms of in- 
visible matter. But it could well be that the discrepancies 
signal a breakdown in our understanding of gravity. 

Among the many theories attempting to extend GR 
some have tried to address the last issue, doing away 
with the need for non- visible or "dark" matter to explain 
anomalies at galactic, cluster and cosmological levels (see 
e-g- (SHI)- Such theories have been labelled "MOdified 
Newtonian Dynamics" (MOND), even though they have 
now been embedded into fully relativistic field theories 
(e.g. of them new effects are triggered be- 

low an acceleration scale, oq ~ 10"^'^ ms~^, a property 
suggested by the phenomenology. The fact that MON- 
Dian behaviour is physically triggered by an acceleration 
scale does not preclude writing ao in terms of length scale: 



io = — 
ao 



(1) 



for which ao could be a proxy. It is interesting that 
Lo ^ 30, 000 Mpc is of the order of the current hori- 
zon/Hubble radius. Nonetheless, it is the onset of low 
acceleration that triggers new effects. MONDian theo- 
ries do not have preferred frames and do not break dif- 
feomorphism invariance; yet new effects emerge in the 



non-relativistic approximation when the total Newtonian 
force per unit mass falls below ao. For this reason one 
may expect the presence of "MONDian habitats" in the 
small regions encasing the saddle points of the gravita- 
tional potential in the Solar System, the points where the 
Newtonian force vanishes [13[. The prospect of a MON- 
Dian saddle test has motivated extensive work fT^-[I^, 
with Lisa Pathfinder (LPF) in mind, but not only. 

One of the most powerful tests of GR, by now elevated 
to the category of "classical" test, employs the echo time- 
delay effect. By flashing "light" (usually a radio wave) 
at a distant object, and catching its reflected "echo" , one 
measures a distinctive delay, if its path intersects a strong 
gravitational fleld. This so-called Shapiro effect was first 
detected with the radar echo off Venus in superior con- 
junction, leading to stringent constraints on the 7 PPN 
parameter |l| . Since then the observational front has im- 
proved very fast. These tests use the fact that for a large 
class of theories (metric theories) the travel time is given 
by 



]dz 



(2) 



where $ is the total gravitational potential, as obtained 
in the non-relativistic limit. Whilst the delay along a 
single path may be gauged away, the variation in de- 
lays along neighbouring paths is operationally meaning- 
ful, and constitutes a bona fide observational target. 

Lunar Laser Ranging (LLR) is a major asset, among 
other fields, in gravitational physics (see for exam- 
ple j20j). Using lunar retrorefiectors one may time very 
accurately echoes of sharp laser signals. It is immediately 
obvious that in principle LLR could probe the Moon sad- 
dle during a Lunar eclipse. Given that the Sun-Earth 
saddle is within the Lunar orbit, it turns out that it can 
also be targeted by LLR during a Solar eclipse. In prac- 
tice a LLR saddle test requires the correct vantage point 
on Earth or even in its orbit. Very Large Baseline In- 
terferometry (VLBI) is another asset in gravitational sci- 
ence (e.g. ^21,5), providing another strategy for probing 
delays across saddles. It correlates images of the same 
object (say, a quasar) as obtained in different continents. 
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A set up could be arranged in which one hght ray goes 
through a saddle whilst the other (a few thousand kilo- 
meters away) does not. VLBI has the advantages that 
it could be used to probe other saddles, (e.g. the rather 
large Sun- Jupiter saddle), and that it relies only on the 
presence of a source behind the saddle, rather than on 
the vagaries of eclipses. In spite of potential difficulties, 
a saddle test with LLR or VLBI could be carried out 
without the logistic overheads associated with the LPF 
saddle extended mission. 

It turns out that the time-delay effect is negligible for 
the MONDian models usually taken as targets for LPF. 
But there are also models well beyond the reach of LPF 
which predict a strong time-delay signal. Therefore the 
experimental test examined in this paper is complemen- 
tary to a LPF test. This is hardly surprising. LPF ac- 
celerometers are sensitive to tidal stresses, i.e. the second 
derivatives of the gravitational potential (they feel their 
variation on a given frequency range). A time-delay test, 
instead, is sensitive to the integral of the potential along 
the line of sight (or rather, to its variation across neigh- 
bouring paths). It is therefore natural to find comple- 
mentarity between the two measurements. An extension 
to LPF would constrain modified gravity theories which 
do away with the need for dark matter. A time-delay 
measurement would instead probe theories which encode 
the property that gravity is subject to confinement, as 
we shall now see. 

It has long been speculated that gravity might resem- 
ble QCD in some aspects. Specifically, gravity could be 
asymptotically free (or "safe") and reciprocally be sub- 
ject to confinement, with a divergent strength at large 
distance, or low energy. The latter is actually what hap- 
pens in the presence of a negative cosmological constant, 
inducing an attractive force satisfying Hook's law (di- 
verging like r). But more generally the issue has been 
raised in the context of the renormalization group flow 
of quantum general relativity (e.g. [2^ - |23 |). where it has 
been conjectured that the degenerate fixed point, when 
lifted, contains divergent IR behaviour. An equivalent 
implementation is possible in the context of MONDian 
theories, should these be released from their obligations 
as dark matter alternatives, but with a key property kept 
and exacerbated. Existing MONDian theories already 
have the peculiarity that they enhance the strength of 
gravity in situations where the standard Newtonian force 
would become weak. Specifically, for astrophysical appli- 
cations, as Fn ^ ao, we have that the MONDian force 
goes like -F^ ~ V^^n, at least in spherically symmetric 
situations. Thus, the MONDian force still drops to zero 
with Fjv, albeit slower. But what if it diverged instead? 
For example, we could imagine the "dual" behaviour 
Fcf, (X 1/F]S[, and once this possibility is considered we 
could consider sharper divergences, such as exponentials. 



with p > very large. If one is to avoid appealing to 
dark matter clearly ao would have to be smaller than the 
usual one, but not otherwise. It turns out that it would 
be very difficult to rule out theories of this sort, except 
for their echo and VLBI saddle delays, as now show. 

We first briefly lay down the formalism for defining 
MONDian theories, without wedding ourselves to a spe- 
ciflc formulation. As explained in [1^], in spite of the 
large number of MONDian theories, the expression for 
the non-relativistic potential invariably satisfles 3 types 
of equations only (which may be formally reduced to 
2). The dynamics may be written as resulting from the 
usual Newtonian potential $Ar and a "flfth force" field, 
(/), responsible for MONDian effects, with total potential 
$ = -\- (j). For "Type I" theories is ruled by a 
non-linear Poisson equation: 



V • {fi{z)V(l)) = kGp, 



(4) 



where, for convenience, we pick the argument of the free 
function function fj, as: 



Ait oo 



(5) 



where k is a dimensionless constant. For "Type 11" the- 
ories we have instead: 

= -J^V • {v{v)V<^n) (6) 
where the argument of free function v is given by 

.^(f)'!^. (7) 

V47r/ ao 

Should these theories serve their duties as dark matter 
alternatives we should require /i ^ z for 2; ^ 1 and 
V ^ I/a/w for w ^ 1, and oq should be the usual MON- 
Dian acceleration. However these theories have an inter- 
est in their own right: they may generally be regarded as 
theories with a preferred acceleration scale. More general 
functions fi or v and values of ag should then be consid- 
ered. If we want to keep the alternative to dark matter 
rationale, then the ao used here must be at least one order 
of magnitude smaller than that employed in traditional 
MONDian theories. However, if we detach these theories 
completely from that role, and if dark matter does exist 
and play a role in the dynamics, this is not true. 

We illustrate our calculations using type II theories, 
because they're simpler. For the purpose of investigating 
confined gravity we shall consider free functions which 
for V 1 are power laws: 



1 



(8) 



F,, 



(3) 



From ([6]) and ([3]) we have p — n — 1 and so for n > 1 
(to be contrasted with the MONDian n = 1/2) we obtain 
confinement behaviour. For n — 3/2 the theory mimics a 
negative cosmological constant in spherically symmetric 
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situations. The dual behaviour suggested in ([3]) with 
p — 1 follows from n = 2. Note that a negative Lambda 
is the perfect dual to standard MOND behaviour. The 
potential (j) diverges around a saddle if n > 2. However, 
this only translates into a divergent time-delay at the 
saddle if n > 3, with ti = 3 representing a logarithmic 
divergence 

Even though parameter k will not appear in the final 
answer for the time-delay, it is important for justifying an 
approximation and setting a scale. The rationale for its 
appearance in ^ and proportionality constant in ^ is 
as follows |14|. If — T' 1 at large argument then Newton's 
constant G is renormalized by k/^Att), and this should be 
small. But so that -F^ ^ ao when Fm ^ ao, we should 
use ([7]) for argument of ly, ii v ^ is to be triggered 

at u 1 in the usual MONDian theory. With (jS]) the 
same requirement becomes: 



|V$ 



ao 



(9) 



with v on 1 /u" triggered at v ^ 1 . The region where field 



goes strongly MONDian is now of size: 



ro 



Oq 

A 



(- 



(10) 



However field (j) is subdominant with respect to the New- 
tonian potential until we get to a distance: 



Oq 

A 



(11) 



from the saddle. It is only inside this inner bubble that 
(j) is both MONDian and dominant. 

For orientation purposes we first evaluate the delay 
effect in Newtonian theory. Introducing cylindrical co- 
ordinates, {z, p, 6}, the delay is obtained by integrating 
Eq. ^ along paths of constant p = b. (This is valid 
for eclipse LLR only; the geometry is more complex for 
VLBI.) In the linear approximation the Newtonian po- 
tential is ^ a(^-^ + For the Earth-Sun saddle 

the tidal stress is A sa 4.6 x 10""s"^. Evaluating ^ 
over z G (— L/2,L/2), retaining only terms that vary 
with impact parameter b, leads to At = — ^ . Even if 
the linear approximation were valid throughout the whole 
flight, we'd get variations of the order of 10~^®s for im- 
pact parameters b ~ 1000km (with a delay at the center 
with respect to outer trajectories). For the Jupiter sad- 
dle the effect is even smaller since A ~ 1.8 x 10~"'^^s~^. 
The Newtonian delay is therefore negligible. 

We now repeat this calculation for type II theories, 
mimicking the calculation of (j) in [l3| , for a free function 
of form ([5]). The potential satisfies ansatz: 



1 



A 



n—1 ^n— 2 



(/0+/2COS(2V^)+/4COS(47A)- 



(12) 



where the parameters / have to be determined numer- 
ically. A particularly simple case follows from n = 4. 



The integration should be performed along z within the 
region of strong MONDian behaviour, delimited by tq as 
defined above. For n — A the integration can be carried 
out explicitly, with leading order term in b/rp: 



At; 



fa . ro 
„ , „■ — arctan — 
c^A^ b b 



(13) 



If we can assume that tq 3> & (always true if k <C 1), this 
becomes asymptotically: 



At: 



(14) 



Parameter k does not appear in the final answer, as long 
as it is small enough to justify taking rg/b — >■ oo. For 
more general n the calculation is more elaborate, leading 
to the asymptotic result {b/ro 0): 



At 



where C is given by: 

^ - - 'n~2\ 1 JO 



C 



hn~3 ^3 



A^- 



(15) 



4 (n-4)(n-6) 
-77/2 H 7 ;:tt J4 



n{n - 



As we see, n = 4 is a particularly simple limit of this 
expression. We note that for n > 3 the relative time vari- 
ation at the saddle diverges, decreasing as a power-law 
in b as we move away from the "bull's eye" . In principle 
the constant C can be positive or negative, leading to a 
delay or an advance at the bull's eye, but we shall call it 
delay for dcfiniteness. 

In spite of this "divergence" the observational implica- 
tions are less dramatic than might be expected. Using 
the cosmological length scale Lq defined in Eqn.(IT]) and 
the strong MOND bubble scale f (in Eqn. lfTT]) ) we can 
rearrange Eq. (jl4p in the suggestive form: 



At« C- 



b f 
cLo 



r 



(16) 



The first factor (besides C) is just the time it takes 
to cross the region closest to the saddle. Unlike time- 
delays caused by the Sun this is small, because the dis- 
tances involved are small: b should be smaller than f and 
f ~ 2.2 m for the Earth-Sun saddle and f ^ 5.5 km for 
the Jupiter saddle. In addition the second factor relates 
the MOND bubble size to the horizon scale, introducing 
a tiny factor. Therefore, even though the third factor 
predicts a divergence, this will happen very close to the 
saddle and be observable only for very steeply diverging 
functions. 

We spell out this expectation in Figs. [T] and [21 which 
describes the situation for the Earth-Sun saddle (as a tar- 
get for LLR during a Solar eclipse) and the Jupiter-Sun 
saddle (as a target for VLBI), respectively. In both of 
these figures we have plotted the (base 10) logarithm of 
the time-delay in picoseconds, as a function of the im- 
pact parameter b (in meters) and the exponent n used in 
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FIG. 1: The LoglO of the delay in picoseconds, as a function 
of impact parameter 6 and exponent n for the Earth-Sun sad- 
dle, as probed, say, by LLR during a Solar eclipse. As we can 
see the delay goes very quickly from very small to very large. 
Realistically, with current technology, only an integrated ef- 
fect might be observable, and even then only for large n. 
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FIG. 2: The LoglO of the delay in picoseconds, as a function 
of impact parameter b and exponent n for the Jupiter-Sun 
saddle, as potentially probed by VLBI. Again the delay goes 
very quickly from very small to very large. Here the bound- 
ary is closer to realistic experimental parameters, and lower 
exponents n come within reach. 



free-function v. In both cases we observe a very abrupt 
transition from the very small to the very large, with 
the contour labeled zero denoting the rough borderline 
for observability with current technology. Typically the 
Earth saddle would have to be probed closer to a meter 
and even then assuming large values of n (in the range 
20-30). The Jupiter saddle might be more forgiving, and 
small values of n ~ 5 could come within reach for b of 
the order of a meter, with n ^ 20 still constrained even 
for impacts of the order of a kilometer. In all fairness 
we cannot be overenthusiastic about the detectability of 
this effect in the first setting, where with current tech- 
nology it would be seen at best as an integrated effect 
(the wavepackets often have a width of about 200 me- 
ters). The next generation of lunar rctroreflectors could 
be necessary. The second situation might be more hope- 
ful. We illustrated our conclusions with the Earth-Sun 
saddle during a Solar eclipse but similar results apply for 
the Moon saddle, as targeted by LLR during a Lunar 
eclipse. Likewise what we have shown for the Jupiter 
saddle has a closely related counterpart with the Saturn 
saddle. Incidentally, Eq. ([T6|) can be used to prove that 
the effect for standard MONDian functions {n = 1/2) 
is negligible {At ~ 10"'^^ s with b about a meter from 
the saddle) . Likewise it can be shown that the functions 
considered here would have negligible effect for a LPF 
test. 



What about other MONDian theories? As an exam- 
ple, we briefly discuss what was labeled type I theories 
in [l4| . For these the situation is more complex due to the 
well known presence of a curl field, softening the diver- 
gence [Tsl , [T^ l . Under strict spherical symmetry this field 
vanishes, so although this is not applicable to a saddle, we 
can gain some intuition. Parametrizing the free-function 
in (g]) as: 



for z ^ 1, we see that ignoring the curl field, the expo- 
nent pin ^ is p = l/(m — 1) , so that m > 1 becomes the 
condition for confining behaviour. Now, m = 3 is equiv- 
alent to a negative cosmological constant A, and m = 2 
leads to a perfect dual {p = 1). However, conclusions 
about the conditions for a divergence at the saddle are 
more subtle, because the magnetic field cannot be ne- 
glected. If this were the case, then 1 < m < 2 would lead 
to a diverging cj), and 1 < m < 3/2 to a diverging de- 
lay. However, now we can only try out the more flexible 
ansatz: 

- -Ci^ifo + /2 cos(2i^) + /4Cos(4^) + ...) (18) 

and search for a solution numerically (using techniques 
presented, in a different context, in j25|). We find for 
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m = 1.1, for example, that a = 1.25 (instead of a = 9, 
expected if the curl field could be neglected). Thus, these 
theories are even more difficult to constrain than type II. 
The situation is similar with type III (which also have a 
curl field). 

It is interesting that something as dramatic as this di- 
vergence can be so elusive. Furthermore, we have only 
solved the problem to linear order, and many questions 
can be raised beyond the scope of the calculation pre- 
sented in this paper. For example, if gravity is confined 
and infrared divergent, as envisaged here, could there be 
a singularity at the saddle? If so, would this singularity 
be naked, or rather, would there be a horizon? Whilst a 
positive answer to the first question is plausible, the an- 
swer to the second question is far from obvious. In both 
cases the detectability of the time-delay effect for the 
free functions used above, as calculated here, is unlikely 
to improve. The field is not attractive, so the usual argu- 
ments about accretion disks and X ray emission do not 
apply (even considering, e.g. the solar winds). One may 
think it odd that naked singularities or horizons could 
be fioating around in the Solar system, but in practice 
the regions where such extreme behaviour is felt are very 
small, and they could pass unnoticed. 

Of course one could fiuff up the divergence region by 
introducing functions of the form: 



(19) 



(v - ly 



in type II theories, for example. Then, the non-linear 
regime would be entered close to the ellipsoid +p^/2 — 
f^, and depending on the details of the full relativistic 
theory, this could signal the formation of a horizon or a 
naked singularity. Either way, assuming LLR geometry, 
any photons with b < 2f would be lost, i.e. they would 
have an infinite time-delay. Close to the disk defined by 
b = 2r the time-delay would diverge as: 



At = C- 



c Ln \ b — 2r 



(20) 



(written in a format to allow easy comparison with ([16])). 
We would now need to be glued to the "horizon" for the 
effect to be measurable, unless n is very large. However, 
we would also have a "black spot", comprising the disk 
b < 2f. This rather extreme free-function is the only 
possibility we found for rendering these theories more 
tangible, and clearly a large cq is then promptly ruled 
out. 

In summary, we hope that in this paper we have 
stressed the radical difference between the gravitational 
physics probed by LLR or VLBI on the one hand, and 
LPF on the other, regarding saddle points. With LPF 
one probes second derivatives of the potential, locally. 
With LLR and VLBI one probes the integral of the po- 
tential, at the end points, as a cumulative effect. There- 
fore with LPF, for standard MONDian functions, we find 
a distinctly changing tidal stress at the saddle (to be con- 
trasted with an essentially DC Newtonian background). 
We cannot realistically get close to the saddle, but even 
far out we can expect signals with large SNRs. With a 
delay test we can potentially probe the region very close 
to the saddle; however, the predicted effects for stan- 
dard MONDian functions are tiny. Nevertheless, we be- 
come sensitive to functions which diverge at low acceler- 
ations, associated with confinement and strong infrared 
behaviour for gravity. Such theories predict extreme be- 
haviour very close to the saddle, raising the possibility of 
singularities. They are beyond the reach LPF and they 
do not purport to present an alternative to dark matter. 
However they have an interest in the own right, and are 
targets for a time-delay test as performed by the current 
or next generation of lunar retroreflectors, and by VLBI. 
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